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1 Introduction 

Exactly soluble models have long provided a foundation for the study of quantum me- 
chanical systems. In the case of quantum systems with one degree of freedom, shape 
invariance (see |]J and the references therein) has proven to be perhaps the most com- 
pelling technique for proving exact solubility and for obtaining the exact solution. In this 
paper, we demonstrate that the Calogero 0] and Calogero- Sutherland |J models — two 
of the best known iV-body exactly soluble quantum models ||. || - exhibit an iV-body 
generalization of shape invariance, and we use this as a tool to analyze some features of 
these models. We will not be able to use this algebraic result to solve completely the 
models in question, but we will outline how this might be done, and make the prelimi- 
nary steps in that direction. Clearly, an algebraic explanation of exact solubility in these 
models would be a significant improvement in our understanding of these models, and so 
our analysis of the shape invariance structure of these models is valuable. 

To lay the appropriate groundwork for our paper, we quickly review shape invariance 
with one degree of freedom in the following section. Then we turn to the heart of the mat- 
ter, obtaining a representation for the Calogero and Calogero-Sutherland Hamiltonians 
in terms of a set of operators similar to raising and lowering operators, and we use this 
representation to demonstrate the iV-body shape invariance of these models. The power 
of shape invariance in the one- dimensional case makes clear the potential significance of 
this result. With this in mind, we then establish a variety of connections, quantitative 
and qualitative, between the Calogero and Calogero-Sutherland models on the one hand, 
and the shape invariant one-dimensional models on the other hand, which provides addi- 
tional evidence that the iV-body shape invariance and one-variable shape invariance are 
indicators of common features in these various models. Our operator representation is 
also valuable as it leads to a natural construction of the supersymmetric extensions of 
the iV-body models we are studying; indeed, we believe this is the first construction of 
the supersymmetric Calogero-Sutherland model. An interesting consequence that we find 
of the iV-body shape invariance is that the supersymmetric extension in each of these 
models is not unique. Note that we examine these supersymmetric extensions, not simply 
because they are interesting in their own right, but because their analysis may well lead 
to the solution of the original non-supersymmetric models, an approach we explore in the 
penultimate section of this paper. 

Readers familiar with shape invariance may be familiar with the result that, for ex- 
ample, the solution of the hydrogen atom may be obtained via shape invariance. Does 
this not represent the application of shape invariance to a problem with three degrees 
of freeedom? The answer is no; shape invariance in this case is actually only applied to 
the effective one-dimensional radial problem. As another example, the multidimensional 
harmonic oscillator completely factorizes into one- dimensional harmonic oscillators, and 
so can be solved by applying shape invariance for a single degree of freedom repeatedly. In 
this paper, we are examining an intrinsically iV-body property, one that does not trivially 
reduce to a one-dimensional problem. 

Recently, Ruhl and Turbiner 11 have explained the integrability of the Calogero and 



Calogero- Sutherland models using Turbiner's definition of exact solubility [0. In this ap- 
proach, given an infinite set of finite dimensional spaces Vo, Vi, V2, ■ ■ ., forming an infinite 
flag Vo C Vi C V% C . . ., an exactly soluble operator h is described as a linear operator 
that preserves the infinite flag of the spaces. With this definition and a series of lengthy 
calculations, they show that the Calogero and Calogero- Sutherland models are exactly 
soluble. Our goal is something different. In the same way that shape invariance has 
illuminated our understanding of exact solubility in one-dimensional quantum mechanics 
(even in those models already known to be exactly soluble), demonstrating that an appar- 
ent analytic coincidence actually arises from an elegant and simple algebraic structure, so 
too do we expect that the algebraic relationships we have uncovered will ultimately lead 
to a ready and efficient determination not only of the exact solubility of these models, 
but indeed of the solutions themselves. Thus we are utlimately interested in the algebraic 
understanding of these models that shape invariance can provide, not simply in the ability 
to solve the models per se. 

2 Shape Invariance for Systems with a Single Degree 
of Freedom 

Why is the condition of shape invariance signficant? In the case of systems with one degree 
of freedom, shape invariance leads not only to exact solubility, but indeed to the exact 
solution. In fact, the algebraic methods of shape invariance probably give the clearest 
understanding of why these models are soluble and certainly give the easiest computational 
methods for obtaining the energy eigenvalues and wavefunctions. We are hopeful that our 
iV-body generalization will be similarly powerful. In order to consider how to proceed in 
the iV-body case, we here summarize the salient features of the one-dimensional case. 

Consider a Hamiltonian H which is a function of some parameter[| a. One wishes to 
find the eigenfunctions and eigenstates of H(a). 

Suppose one can write H(a) = A\a)A(a). As is well known from studies of supersym- 
metry, the partner Hamiltonian H(a) = A(a)A\a) has essentially the same spectrum. 
To be more precise, if 

Hi)) = Ei/j , 



then 

and likewise, if 

then 



H(Ai/j) = A(Hip) = E(Atp) 
H(A^) = E{A^) . 



Thus, the spectrum of non-zero energies is the same for the two Hamiltonians, and the 
corresponding states that have the same energy are easily constructed from each other by 



1 a may as well stand as a collective symbol to indicate dependence of H on many parameters. 



acting with A or A\ The only place they may differ in their spectra is that, for example, 
H may have a zero-energy eigenstate while H does not. Note too that the representation 
in terms of A and A' ensures that the two Hamiltonians are positive semi-definite, which 
means a zero energy state, if it exists, is the ground state. Furthermore, such a zero energy 
eigenstate is especially easy to construct, as one need only solve Aip = 0, a first-order 
rather than a second-order differential equation. 

Shape invariance exists when the two partner Hamiltonians describe the same physical 
model at different values of the coupling constant (s), i.e., when they are related in the 
following way, 

H(a)=H(a 1 ) + R(a 1 ) , (1) 

where a.\ = f(a) and R(a) are both a-dependent constants. When this occurs, one can 
obtain the entire spectrum of H and all the wavefunctions once one knows the ground 
state. Here is how it works. 

The ground state of H(a) we find by A(a)ipo(a) = 0. This has zero energy, and hence 
no "tilde" partner. The first excited state of H is degenerate with the ground state of H. 
However, H{a) = H(ai) + R(ai), and so its ground state is ipo((*i), with energy equal 
to Ei = R(ai). Then the first excited state of H(a) is A' (a)ij)o(&i) an d also has energy 
E 1 = R(ai). Iterating this process, one obtains the exact energies 

( , if k = , 

Ek = \ E R( aj ) , if k > . ( 2 ) 

and corresponding energy eigenfunctions 

ip n (x;a ) = A*(a )A*(ai)...A*(a n - 1 )ipo(x;a n ) , (3) 

where otj+i = f(cxj). 

An example of a shape invariant model is given by the trigonometric Rosen-Morse 
model, with Hamiltonian 

H = p' + ^ .-*. (4) 

sin (ax) 

The parameters of this potential are a = (b, a). Define 

A = — b cot(ax) , (5) 

ax 

A f = — b cot(ax) . (6) 

ax 

Then A' A = H(b, a) and AA^ = H(bi, Oi) + R(bi, ai) where b\ = b + a, a\ — a and 

R(b 1 ,a 1 ) = (b + af-b 2 . (7) 

The spectrum of the model is thus given by 

E n = [b + naf-b 2 , (8) 

where n is a natural number n — 0, 1, . . . . It is worth noting that at b = a, the Hamiltonian 
H is that of a free particle in a box. 



3 TV-Body Models: Exact Solubility and Shape In- 
variance 

The use of shape invariance, or some suitable generalization therof, to solve systems with 
multiples degrees of freedom remains an unrealized goal. Where shape invariance has been 
applied to such systems, for example the 3-dimensional hydrogen atom ||, it has actually 
been applied to the effective one-dimensional radial problem that arises in a sector of fixed 
angular momentum. 

Nonetheless, the idea that there should be an algebraic explanation for the exact 
solubility of systems with multiple degrees of freedom remains compelling. Indeed, it 
seems unlikely, given our experience with ordinary shape invariance, that exact solubility 
should arise simply as an analytic accident. 

Two of the most studied iV-body integrable models are the Calogero and Calogero- 
Sutherland models. In this section, we demonstrate that these models are shape invariant, 
or, to be more precise, satisfy an algebraic identity that is the natural iV-body general- 
ization of shape invariance. We expect that the algebraic results we have found will lead, 
ultimately, to a ready derivation of the integrability of these models, and to the efficient 
algebraic calculation of the actual wavefunctions and energy levels of these systems. What 
we will do in this paper is to obtain these identities, and then discuss the analysis of these 
models that our representation permits. 

The generalization of shape invariance we obtain is as follows. In the models in 
question, we find we can define a set of operators A4 and A] such that the Hamiltonian 
can be written as 

H(a) = A\Ai + ■■■ + A ] N A N . (9) 

However, the operators A4 and Aj do not commute. Thus, this identity provides an elegant 
representation for the Hamiltonian, but nonetheless one that generalizes, but falls short 
of, separability. 

The compelling result, however, is the following algebraic structure which turns out to 
be present in these models. If we interchange the Aj's and Aj's, we obtain the associated 
Hamiltonian 

H(a) = A X A\ + ■■■ + A N A ] N . (10) 

We are able to identify these models as shape invariant because, it turns out that 

H(a) = H(a 1 ) + R(a 1 ) . (11) 

This identity generalizes the notion of shape invariance familiar from models with one 
degree of freedom. We conjecture that the ability to write the Hamiltonian in these two 
ways underlies the exact solublility of these models. 

In the rest of this section, we will present the specific forms of these representations. 
Then, in this and later sections, we will discuss the applications of this representation 
of the Hamiltonians, providing support for our conjecture, keeping an eye toward how 
this iV-body shape invariance should be related to the exact solubility of the models in 



question. Also, in appendix |A|, we present a general definition and treatment of "shape 
invariant" potentials for iV-body potentials in an analytical framework. 

3.1 The Calogero Model 

The Calogero model is an iV-body non-relativistic quantum mechanical model given by 
the Hamiltonian 

ffc—Eaf + EE'r-^iS' (12) 

i ij ( x i~ x j) 

where the prime in the summation sumbol means 'restricted sum', i.e. summation over 
the dummy variable j excluding the value i for which the denominator vanishes. We can 
obtain a useful representation of this Hamiltonian as follows. 
Let us define the prepotential^ 

^ = £'-^- , (13) 

in terms of which we define the operators 

Ai = di + Wi, A\ = -d t + Wi . (14) 

These operators obey the following algebra: 

[At.Aj] = 0, [4,AJ] =0, (15) 

Al.A,] = -2 Si ^=( El 2 l^F' H = ]• (16) 



Two further relevant identities are that 

i 

and that the curvature associated with the prepotential vanishes, i.e., 

diW^djWi . (17) 

The first important observation is that the Hamiltonian can be written in a simple 
way using the Aj's and A|'s. In particular, one finds that 



i? c («) = E4(«)A 



a] 



The quantity we term the prepotential is sometimes referred to as the "superpotential" in the lit- 
erature. We prefer to use the term "prepotential," reserving "superpotential" for the quantity whose 
gradient gives Wi, in order to be consistent with the terminology of supersymmetric field theory. 



at coupling constant 

g = 2a (a — 1) . 

In order to verify this result, note that 

A\A i = -d 2 i +Wf-d i W i . (19) 

To see that this expression yields the Calogero Hamiltonian, it is necessary to verify that 
the three-body interactions contained in the W 2 term cancel. (No other term contributes 
to a direct three-body potential.) To see this, we consider all terms in A\Ai that involve 
Xi, Xj, and Xk, where these three coordinates are distinct. The terms are 

1 1 1 

+ 1 w T + 



yXi Xj )\X-i Xk} \^j *"i)v"j X'k) \Xk Xi)\Xk ^j) 

which it is easy to verify is zero by putting everything over a common denominator. 

This representation for the Hamiltonian in terms of A48 and A- 's is potentially rather 
powerful. Unfortunately, these operators do not diagonalize the Hamiltonian, due to 
the non-commutation contained in (|I6D. Nonetheless, the vanishing of the three-body 
interactions, the representation of the Hamiltonian as a sum of quadratic terms, and the 
pairwise commutation of the Aj's (respectively, the Aj's) does at least give us a structure 
akin to separation of variables. It is the rare Hamiltonian that takes this form. 

There is however a more significant algebraic identity at work here. If we switch A{ 



and A], we obtain 



In other words, 



H c (a) = J2 A t {a)A\{a) = H c (a - 1) . (20) 



£ A(a)Al(a) = J2 4(« + l)A-(« + 1) • (21) 



Interchange of the Aj's and A]'s has the same effect as shifting the coupling constant! This 
is the hallmark of shape invariance, and in the later sections of this paper, we consider 
the possible connection of this identity to the exact solubility of the model. We have thus 
have demonstrated that the Calogero model exhibits an iV-body generalization of shape 
invariance, the first time it has been found in such a model. 

These operators can be used to represent physically relevant quantities. An elementary 
but useful result is that the total momentum can be written as 

iW = -*E^ = *E4, (22) 

i i 

since J2i W* = 0. The quantity Ptot commutes with each of the A4 and A\ individually, 
and hence we can classify solutions in terms of their total momentum (not surprisingly) 
in a way that does not interfere with our operator representation of the Hamiltonian. 

The Calogero model has an important variation which we now describe. One can in- 
troduce additional pairwise interactions among the particles quadratic in the coordinates: 

^ = -E^EE't-^.+EE'^ 2 (^^ + - (23) 



In order to distinguish this from the model we discussed initially, we shall coin the term 
"harmonic- Calogero model" to describe the theory with this quadratic interaction. The 
discussion presented above in the case of the ordinary Calogero model can be repeated 
nearly verbatim, with only minor changes. For simplicity, we here highlight only the 
important details. 

The operators we must define are 



a 



M = $-£'^^+0 £'(*«-*,), (24) 

3 %t X 3 3 

4 = -^-E'^r + ^E'^-^'), (25) 

where 



■ Xi Xj 



2a(a-l)=g , (3 = -^= . (26) 

2x/N 



With these operators, we have both the generalizations of separation of variables and of 
shape invariance that we found in the ordinary Calogero model. To be precise, 



Hhc(a) = E 4(«)A(«) , (27) 



and the constant c in equation ( p^) is 

CO 



*N(N-l)(aN+l). (28) 

Then 

J2 M<*)4(<x) = E 4(«0A(«i) + R(ai) , (29) 

i i 

where 

ai = a+l, R(ai) = -^= VN (N - 1) («i - a) N . (30) 

v2 

Further analysis of this model can be performed exactly as for the ordinary Calogero 

model; we leave the particular calculations to the interested reader. 

3.2 The Calogero-Sutherland Model 

Our formulation of the Calogero-Sutherland model, and the operators we use to rep- 
resent the Hamiltonian, can be developed very similarly to the Calogero model. The 
Calogero-Sutherland model is given by an iV-body, non-relativistic quantum mechanical 
Hamiltonian, 

tf^H-E^ + EE'-Trf^ — vf + c - ( 31 ) 

i i j Sm HXi-Xj)] 

In the following, we shall set a = 1. We now define a prepotential 

W t = -a J2 ' cot(x< - Xj) . (32) 



This leads us to introduce the operators 

A i = d i + W i , 4 = -di + Wi . 



(33) 



These operators satisfy a variety of algebraic relations, and are useful in analyzing the 
Calogero-Sutherland model. For these new A4 and A] operators, we have 



|Aj, Aj 
A-i , Aj 



A i , Aj 



-2diWj 



Sfc 



-2a 



if i = j, 



sm 2 (xi-x k ) ' 



(34) 
(35) 



sin (xj— Xj) ' 

It is worth noting again here, as in the previous model, that 



and that 

Setting 

and 



c = —a 
the Hamiltonian can be written as 



diWj = djWi . 



g — 2 a (a — 1) 

2 N(N* - 1) 



ffcs(a)=EA t («)A 



a 



(36) 



(37) 



Just as with the Calogero model, the three-body interactions automatically cancel. As 
before, the only source of three-body interactions is the W 2 contribution to A\Ai. Let us 
consider all the terms that include Xi, Xj, and x^. For simplicity, let a, b, and c represent 
the differences Xi — Xj, Xj — Xk, and Xk — Xi, respectively. Then we have terms 

— cot a cot b — cot b cot c — cot c cot a . 

Now, since a + b + c = 0, we see 

cot a cot b + cot a cot c + cot b cot c = 1 . (38) 

Thus we can write the Hamiltonian in terms of the A^s and A^s in a form that generalizes 
separation of variables. 

And, again as before, there is more. If we switch the hermitian conjugates, we get the 
partner Hamiltonian 

Hcs(a) = J2 Ma)Al(a) . (39) 



Remarkably, as with the previous model, we see that interchanging the hermitian conju- 
gates is equivalent to shifting the coupling constant! To be precise, we find that 

H cs (a) = HcsM + R( ai ) , (40) 

with 

c*i=a + l, R(a + 1) = (a\-a 2 ) ~ • (41) 

o 

Put it another way, 

]T A(a)A\(a) = J2 A(<*i)M<xi) + R(<*i) ■ (42) 

i i 

We have therefore obtained a generalized shape invariance in another iV-body model. 
Note, too, that we again find 

Ptot = -i ]T A, = +i^ A > 

i i 

where of course Ptot commutes with all the individual A^s and Aj's.. 

Notice that in both these models, the key to shape invariance is that, up to an additive 
constant, diWi and WiWi have the same functional dependence, which is the same as that 
of the potential. This relationship, we expect, and the generalized shape invariance it 
implies, truly underlie the solubility of the system. 

4 Generalizations of 1-body Shape Invariance 

The Calogero and Calogero-Sutherland models share a number of features in common with 
familiar shape invariant models of one degree of freedom. These similarities go beyond 
the algebraic identities above, and relate to the features of the model. We review here 
some of these similarities, laying the groundowrk for a future solution of these models by 
means of shape invariance. 

1. Interchanging Aj and A\ has the same effect on the Hamiltonian as shifting the 
coupling constant. This is the fundamental hallmark of shape invariance in the one- 
body problem, and arises as well in the iV-body Calogero and Calogero-Sutherland 
models. 

2. These iV-body models are known to be exactly soluble. Likewise, the shape invariant 
models with one degree of freedom are known to be exactly soluble (indeed, were 
generally known to be so before shape invariance was discovered). 

3. The ./V-body models in question are the obvious multiparticle generalizations of 
known shape invariant models with one degree of freedom. More precisely, the iV- 
body models have potentials which involve only pairwise interactions; the pairwise 
potential in each case is given in terms of a single relative coordinate; and this two- 
body potential is just one of the shape invariant, exactly soluble potentials from 
ordinary quantum mechanics with one degree of freedom. 

10 



4. The partition function computed semiclassically is exact in these theories. The 
evidence for that has been presented in . This is also a feature exhibited by shape 
invariant models with one degree of freedom 0. 

5. For the case N = 2, i.e., the two-body problem, we actually directly recover the 
conventional case of shape invariance with one degree of freedom, once we separate 
out the total momentum Ptot- 

6. The ground state is easy to construct, by requiring it to be a state annihilated by 
all the Aj's. This gives us a zero energy state, but not for the partner Hamiltonian. 
Shape invariance, however lets us get the ground state for the partner. 

It is worth turning our attention to the last two of these items in particular. We turn 
first to construction of the ground state. 

To find the Calogero or Calogero-Sutherland ground state, we use the representation 
of the Hamiltonian as H = Yi A-^-i- Since this expression is positive semi-definite, the 
states of the theory must have non- negative energy. Therefore, if we can find a zero-energy 
state, this will of necessity be the ground state. 

In fact, since the Hamiltonian is the sum of iV postive semi-definite terms, in order 
for a state $o to have zero energy, it must satisfy 

A&o = , V i . 

It is easy enough to solve this, obtaining 

$o(a; 1 ,...,x n ) = 0oe-^r^^. 

As long as this function is normalizable, it yields a sensible state, and then this is the 
ground state wavefunction of this system. 

Notice that normalizability forces us to recognize that if there are solutions to the 
equation Ai(a)Q (a) = 0, there are no corresponding solutions &o(a) to Al(a)<& (a) = 0. 
Nonetheless, we can find the ground state $'(a) of Yi Ai(a)A\(a) by exploiting the shape 
invariance of the model. In particular, since Yi ^-i( a )A( a ) = YiA( a + l)Ai(a + 1) + 
R(a + 1), the ground state $'(«) oc $(a + 1). The energy is given by the constant 
R(a + 1). Thus shape invariance allows us immediately to obtain the ground state energy 
of Yi AiA], despite the absence of a physical solution to AJ$ = 0. 

The second point we wish to turn to is the special situation when N = 2, i.e., when 
we have a 2-body Hamiltonian. In our two models in this case, then, the Hamiltonian is 

H^ = A\A X + A\A 2 , (43) 

where the superscript indicates the number of particles. We can write the Hamiltonian 
in this case as 

H {2) = \{A + A)(M + M) + \{A - A){M - A 2 ) . (44) 

11 



Recall that the total momentum 

Ptot = -i(Ai + A 2 ) = i(A\ + 4) (45) 

and that both A\ — A 2 and A\ — A\ commute with Ptot- Consequently, if we define (for 
simplicty of notation) B = Ai — A 2 , we can write the Hamiltonian as 

# (2) = \p 2 tot + \&B . (46) 

Since [Ptot, B] = [Ptot, B^] = 0, the solutions of this model are 

^{x l) x 2 ) = e lk ^ +x * ) ^ B {x l -x 2 ) , 

where iPb(xi — x 2 ) is an eigenstate of B^B. 
What of the partner Hamiltonian, 

H 2 = A 1 A\+A 2 At ? (47) 

Using the same definitions as above, we see that 

#2 = \Ptot + BB^ (48) 

Thus, other than the overall center-of-momentum term exp[ik(xi + x 2 )\ in the wavefunc- 
tion, solving this problem simply involves finding the wavefunctions of B^B and BBK 
This can be done using the familiar result for shape invariance for systems with one de- 
gree of freedom, and thus we see that our identity in the iV-body case does indeed reduce 
to the familiar shape invariance, with all its attendant usefulness for exact solubility, in 
the two-body case. Later in the paper we explore some ideas of how our identity might 
produce an exact solution in the iV-body case. 

To close this section, we note that there is a generalization of our separation in the two- 
body case into Ptot and B. While unfortunately this generalization does not lead directly 
to the exact solution of the problem, it does give a representation of the Hamiltonian in 
an alternative basis that is superior in certain respects to the representation in terms of 
the Aj's and A^s above. 

Let us go to a Jacobi 3 basis. We define 

Bi = ^(A 1 -A 2 ), (54) 



3 Our terminology comes from the similarity of the above transformation to the choice of Jacobi 
coordinates. In classical mechanics, the Jacobi coordinates 

(49) 
(50) 

(51) 
x N , (52) 

(53) 
mi +m 2 + ■ ■ ■ + rriN 

are used to separate out the center-of-mass motion. 

12 



2/1 = 


= xi - x 2 , 


2/2 = 


mill +m 2 x 2 
mi + m 2 




mi x\ + m 2 x 2 + ■ ■ ■ + m N -i Xn-i 


Vn—i 


mi + m 2 + ■ ■ ■ + m N -i 


Y = 


mi x\ + m 2 x 2 H + mn xn 



B 2 = -L(A 1 + A 2 -2A 3 ), (55) 

(56) 

B N ^ = , [A 1 + A 2 + --- + A N _, - (TV - l)A N ] , (57) 
y/N(N - 1) 

Sjv = -^(Ai + Aa + .-. + Ajv). (58) 



In this operator basis, we still have the generalized separability result, 

H = Y / BJB l 

i 

and we can still see directly the TV-body shape invariance result 

J2 B l BJ(a) = J2 BJB^a,) + R( ai ) . 

i i 

The advantage of this representation is that the commutation relations are, though less 
symmetric, in one respect, a little simpler. We again have 

[B i ,B j ] = [BlB}]=0, 

for all i and j. But now the operators B^ and B' N are proportional to the total momentum, 
and so we have 

[B N ,B}] = [B* N ,B j ] = 0, 

for all j. In addition, this choice automatically pulls out the total momentum, while 
preserving a raising/lowering operator formalism. In other words, we can write 

1 N ~ 1 

H = 2 P TOT + Yl B i B i 

~ i=l 



and 



___ 1 N ~ l 

H = -PIot + E BiBt 
z i=i 



5 Super symmetric Construction and Generalizations 

For the case of shape invariance with one degree of freedom, the algebraic results are best 
understood in the context of supersymmetry. Thus we anticipate that the supersymmetric 
generalizations of the Calogero and Calogero-Sutherland models offer the best possiblity 
for understanding the exact solubility of these models. To this end, we here review the 
relevance of supersymmetry to shape invariance with one degree of freedom, and then 
proceed to construct the supersymmetric Calogero and Calogero-Sutherland models. We 
will see that our representation of the Hamiltonians in these cases leads to an easy and 

13 



natural construction of supersymmetric generalizations for these models. In fact, we be- 
lieve that this is the first appearance of the supersymmetric Calogero-Sutherland model in 
the literature; the supersymmetric Calogero model has been discussed in |§ [UJ. Inciden- 
tally, we mention that recently higher dimensional supersymmetric quantum mechanics 
has been studied by other authors too |ll[, although with a different perspective. 



Let us first therefore consider quantum mechanics with one degree of freedom. The 

degeneracy of the spectra of the two Hamiltonians H and H can be understood as a 

consequence of supersymmetry. One constructs a single supersymmetric system, with 

Hamiltonian 

H 0/ 

_0 H_ 

Defining the operators 



n 



(59) 



we can easily verify the algebra 



Q = 


'0 0" 
A 0. 


•> 


Q f = 


■o At 

.0 


•> 


a. 

[n,Q] = 


[n,Q^} = o, 


{Q,Q} = 


{QW} = o 


{Q,Q ] } = 


n . 





(60) 
(61) 



(62) 
(63) 
(64) 

The ground state, if it has zero energy, is a supersymmetry singlet; other than that all 
states must be paired. 

Recognizing that a^ measures fermion number in this example, one sees that using 
shape invariance and supersymmetry to link together wavefunctions at different fermion 
number is what enables one to solve the theory. Motivated by this one-body result, we 
now construct the supersymmetric generalization of the iV-body Calogero and Calogero- 
Sutherland models, and then study the ways in which the wavefunctions and energy levels 
at different fermion number are related to each other. We will find that going from fermion 
number to the maximal fermion number N, one travels from the original Calogero or 
Calogero-Sutherland Hamiltonian to its shape invariance partner. While we will not be 
able to relate these two fermion number levels directly, we will identify an approach for 
obtaining such a relationship, and hence solving the original models via shape invariance. 

5.1 Supersymmetric Generalizations 

One advantageous feature of our representation of the Hamiltonian and related physical 
quantities in terms of A^s and Aj's is that it makes construction of the supersymmetric 
generalizations of the Calogero and Calogero-Sutherland models completely straightfor- 
ward. We thus can and do use the same formalism for both models. 
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Let us define Grassmann variables ipi and ipl with the anticommutation relations 

{^i,^} = {4,^} = 0, (65) 

{ipi , i)]} = 5y . (66) 

We define the supercharges 

i i 

and the Hamiltonian 

It is easy to verify that Q 2 = {Q^) 2 = 0. Thus this theory has aniV = 2 supersymmetry; 
both the real and imaginary parts of Q square to the Hamiltonian. 

The structure of these super symmetric theories is related to the generalized shape 
invariance. We note that the superalgebra determines that all states are annihilated by 
Q or Q^ . The zero energy states are annihilated by both. Of the remaining states, those 
eigenstates of the Hamiltonian annihilated by Q are also in the range of Q; and those 
annihilated by Q^ are in the range of Q^ . 

5.2 Relating States at Different Fermion Number 

Ignoring the zero energy states for simplicity, we consider dividing the space of eigenstates 
into states of definite fermion number (which can be anything from to N) and also 
according to whether they reside in ker Q or ker Q\ 

1. The 0- fermion states: These states are automatically in ker Q. Furthermore, they 
satisfy the equation 

Y,A l {a)A\{a)<P = E<p. (67) 

i 

In other words, the states of the ordinary bosonic Calogero or Calogero- Sutherland 
model are the 0-fermion states of the supersymmetric version. 

2. The iV-fermion states: These states are automatically in ker Q\ Futhermore, they 
satisfy the equation 

Al(a)Ai(a)(f) = E<\> . (68) 

In other words, the iV-fermion states are also the states of the ordinary bosonic 
Calogero or Calogero-Sutherland model — except at a shifted value of the coupling 
constant, due to shape invariance. 

3. The 1-fermion states in ker Q^: These are obtained by applying Q< to the zero 
fermion states. Thus, these states are degenerate with and derivable from the 0- 
fermion states, i.e., from the original Calogero or Calogero-Sutherland model. 
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4. The 1-fermion states in ker Q: The best way to understand these is as follows. Let 
us write them as 

J2 &(zi,---,aJjvM|0) • 

i 

Then it is straightforward to show that A^i = 0, and that 

J2 (AAUj - A t A](f) t ) = \(f) j . 

i 

Let us define $ = J2i4>i- One immediately determines that since J2iA oc Ptot 
commutes with the Ai, 

AiAl$ = \$ . 

This means that $ = J2i 'Pi is either a 0-fermion wavefunction (i.e., a solution of 
the original Calogero or Calogero-Sutherland model) or is itself vanishing. Thus, as 
long as the sum does not vanish, the energy of such a wavefunction is degenerate 
with a state of the zero-fermion spectrum. 

5. The 2-fermion states: Those in ker Q* we obtain from the 1-fermion states in ker Q. 
The others we construct similarly to what we did above. 

6. (N — l)-fermion states in ker Q: These are simply obtained by applying Q to the 
iV-fermion states, and so we get a spectrum and set of eigenf unctions directly from 
the ordinary Calogero or Calogero-Sutherland model at shifted coupling constant. 

7. (N — l)-fermion states in ker Q^: We can proceed much as we did for the 1-fermion 
states in ker Q. We find that the (N — l)-fermion states in ker Q^ can be written 
as 

i 

where | N ) is the state with iV fermions in the fermionic subspace (i.e., proportional 
to ifti ■ • ■ ipN | )), and x( x ii ■ ■ ■ > x n) is the bosonic piece of the wavefunction. Then 

Axi = o 

and one finds too that X = J2i Xi either vanishes or is itself an eigenstate of J2i A]Ai, 
and so the energy spectrum includes the energy eigenvalues of the bosonic Calogero 
or Calogero-Sutherland model at shifted coupling constant. 

5.3 Comments 

We see from this method that we can related the spectra at different fermion number. 
Ultimately, we expect that refinements in this procedure will enable one to relate the 
spectra of the 0-fermion and iV-fermion states, which means relating the Hamiltonians H 
and H, which in turn means relating the Hamiltonians H(a) and H{a,\). Consequently, 
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we would expect to be able to solve the model algebraically, just as one does in shape 
invariant systems with one degree of freedom. 

A brief look at the case N = 2 is worthwhile. While we discussed the case N = 2 
without supersymmetry and saw how to use generalized shape invariance to solve such a 
system, we see here how supersymmetry organizes the states usefully. In the case N = 2, 
the 1-fermion states in ker Q for example are degenerate with and can be obtained by 
acting on the 2-fermion states with Q; but they may also be related to the 0-fermion 
states as described above. The bosonic part of the state has components <\>\ and <p2, and 
as long as these do not cancel additively, the sum <p\ + 02 is a 0-fermion state that is 
degenerate with the 1-fermion state, and hence in turn with a 2-fermion state. 

The case of iV = 3 also deserves special study, as everything is still very tightly 
constrained. Here, the spectra of the 1-fermion states in ker Q and the 2-fermion states 
in ker Q^ are degenerate; indeed, using the labels 4>i an d Xi as above, one finds 

<pi = ±c^2 tijkAjXk (69) 

and a conjugate expression for \i- Since these spectra are degenerate, and in turn are 
related, respectively, to H(a) and H(a,\) as discussed in items 4 and 7 in the preceding 
subsection, one has nearly all the ingredients in place to relate directly the spectra of 
H(a) and H{a,\) in the N = 3 case. 

5.4 Non-Uniqueness of Supersymmetric Extensions 

Exactly because the Calogero and Calogero-Sutherland models are shape invariant, their 
supersymmetric extensions are not unique. This is easy to see, and important to recognize. 
Let us consider the pure Hamiltonian 

H=-8% + gV + c, 

whereQ g = la (a — 1). We can write this Hamiltonian in two distinct ways, 

H = J2Al(a)A i (a) 

i 

and 

# = E A t (a + l)A\(a + l) . 

i 

Up to an additive constant, these Hamiltonians are the same, and thus describe exactly 
the same physics. We can take either one of these representations as the starting point 
for our supersymmetric construction. However, although both these models will have the 
same 0-fermion sector, they will have different states at higher fermion number^. 

4 In this subsection, we have separated the coupling constant dependence from the potential. 

5 This is in addition to the more obvious observation that one can consider two generalizations of the 
original Hamiltonian, one in which it appears in the zero-fermion sector and one in which it appear in 
the ./V-fermion sector. This, however, is unrelated to shape invariance, and merely involves switching ipi 
and ipj. 
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Thus, because of the iV-body shape invariance, we have the two supersymmetric 
Hamiltonians, which have the same bosonic sector but are otherwise distinct, 

Hsi = {At(a)4,M^i}, (70) 

and 

H S2 = {A^a + 1)V4, 4(« + 1)^} • (71) 

6 Ending Remarks 

We have identified the algebraic identity that we expect underlies the exact solubility of 
the Calogero and Calogero-Sutherland models. This identity is a natural iV-body general- 
ization of shape invariance, and we have explored its implications for the analysis of these 
models, indicating in particular how it relates to integrability. We have also seen how 
this identity can be incorporated directly into our understanding of the supersymmetric 
Calogero and Calogero-Sutherland models; indeed, it appears that the supersymmetric 
context will be the natural arena for using the generalized shape invaraince to solve the 
original non-supersymmetric models. Thus, the work in this paper should provide a foun- 
dation for obtaining a simpler understanding of the solutions of Calogero and Calogero- 
Sutherland models, and a deeper understanding more generally of exact solubility in the 
TV-body case. One particular point to mention is that the eigenfunctions of the Calogero- 
Sutherland model are deeply related to the so-called Jack polynomials of the mathematics 
literature. These polynomials are defined in the space of symmetric polynomials obeying 
certain conditions ]ll| • Although, they are quite well studied, mathematicians are lacking 
an explicit representation. Recently, Lapointe and Vinet presented such a formula |p~3|| . 
Our proposal for an extended shape invariance principle in higher than one dimensions 
might prove useful to either derive this formula in a compact and natural way or find an 
alternative simpler one. 
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A Theory of Many-Body Shape Invariant Potentials 

Here, we systematize the theory of iV-body shape invariant potentials. We develop the 
theory in a completely analogous way to the 1-dimensional theory so that the similarities 
and differencies are apparent. 

A.l Factorization of the Schrodinger Equation in the Many- 
Body Case 

Let us consider a potential V (xi, x 2 , ■ ■ ■ , x^) with ground state wavefunction \Pq (x\, . . . , xjy) 
and ground state energy Eq . As in the 2-particle case, without loss of generality, we can 
assume that Eq — 0. Now from the time-independent Schrodinger equation 

H<® ^ = E^ ^ 0) = , (72) 

where 

N 

if(°) = -Y / d*V (x 1 ,x 2 ,...,x N ) , (73) 



i=l 

where <9« = d/dxi. 

Obviously, the potential Vq{x\, x 2 , ■ ■ ■ , x^) is related to the ground state wavefunction 
^o by the equation: 

N -)2 lT ,(0) 

(74) 
mihilation operators: 

(75) 
(76) 

(77) 
(78) 



V (x) = 


N 02^(0) 
i=l ^0 


In this case we introduce, the following 


creation and 


A = 


f) i ° 

1 ^ 


A t = 


d 9^ 


We clearly see 

H® = 


N 


Finally, in the many-body case we define 


Wi = 


o^P 



In terms of the functions Wi, the creation and annihilation operators take the form 

A\ = -di + Wi, (79) 

A = d l + W i . (80) 
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A. 2 A Special Class of Potentials 

A special class of potentials that has been discussed thoroughly in the literature is the 
class of potentials for which the ground state can be written in the form 

N N 

< ] = n n 4 0) (^-*i)- (si) 

i=l j=i-\-l 

This demand leads to the condition 

A djjj^jxi -Xj ) 

Pi l/>o{Xi-Xj) 

If W(x) is the corresponding prepotential for the 2-body problem, we see that 



Wt - E (o°) -T ■ ( 82 ) 



N 

Wi = J2' W ( x i~ x j) ■ (83) 

3=1 

Now, the previous choice of the wave function determines uniquely the potential 
Vq{x\, x 2 , ■ ■ ■ , Xpf). Using relation (|74|) , we find that 

ir , x v^ dft/j (xi -Xj) ^ diif) (xi - Xj) diip {xi - %k) 
V G {x l ,x 2 ,...,x N ) = l^—- ( ^-+ 2^ 



j# ^o(Xi - Xj) i^k^i ^(Xi - Xj) 4> (xi - x k ) 

= 2\2 v o( x i - x j) + 2\1 W ( x i ~ x j) W ( x i ~ x k) ■ (84) 

We notice that the last term consists of a sum of terms 

W(a - b)W(a - c) + W(b - a)W(b - c) + W(c - a)W(c - b) 

Let A = a — b, B = b — c, and C = c — a for which A + B + C = 0. Although we are 
not restricting ourselves to shape invariant potentials at this point, we make the following 
observation. For the usual shape invariant potentials, except the Morse potential, W{x) 
is an odd function; so it is quite useful to restrict ourselves to odd super-potentials. 
Therefore, the expression written above has the form 

-W(A)W(C) - W(A)W(B) - W(C)W(B) . 

One now sees that direct 3-body interactions are avoided if there is a function Vq{x) such 
that 

- W(A)W(C) - W(A)W(B) - W(C)W(B) = v (A) + v (B) + v (C) . (85) 

In this case the potential Vq is written 

Vq(xx, ..., x N ) = 2\2 v o( x i ~ x j) + 2\2 vo(xi - Xj) . (86) 



There are a few known solutions [14]] to equation (]85|): 
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W(x) 


v {x) 


Vo(x) 


^o(x) 


X 


b(6 t 1} +a 2 x 2 + 2ba + a 


ab + yi 2 


|™|6 „ax' 2 /2 


asgnx 


a 2 sgn 2 x + aS(x) 


a* 

3 


g o|x| 


a cot a; 


a(a-l) 2 
■ 2 W 

sm x 


3 


sinx a 


a coth x 


a(a-l) . o 

. , 2 -f- CI 
sinh x 


a 2 
3 


sinhx a 


a>C( x ) ~~ ^^- x 


W 2 (x) - W'{x) 


l/2P(x)-iy 2 (a;) 


|^i(7rx/2a;|ir/2a;)| 



In the above table, £(ar) is the Weierstrass zeta function, C'( x ) = ~'P{ x )-> an d oj is the 
half-period (15 . 



A. 3 Shape Invariant Many-Body Potentials 

We now introduce the depence of the potential (|86|) on possible parameters a : 



N N 

J2 M x i - x v a o) + J2 v ( x i - Xj] a ) , (87) 



Vo(a;i,...,xjv;ao) 

As we have discussed, the iV-body hamiltonian is written as 

TV 

# (0) ("o)=E4(«o)^(ao) ■ 



We notice now that the operators A(a) and A\a) satisfy the commutation relations 

[Ai(a), Aj(a')} = +W'(xi - xf, a) - W'(xj - Xi\ a') , 
A\(a), A](a')] = -W\xi - xf, a) + W'(xj - x { ; a') , 

Aj(a:), A^ct')! = —W'(xi - xfa) —W'(xj - Xi\a') , 
if i 7^ J and 



[A(aO,A(aO] = -J2' w '( x i- x ^ a ) + Y,' W '( x i- x ^ a ') > 

fc fe 

4(a), 4(a')] = +E' W/ '(^-^;a)-E'W / '(^-^;«') . 

fe it 

A(a),4(a')] = +E'^'(x J -x fc ;«) + E , ^ / (x J -x fc ;a / ) • 



(89) 
(90) 

(91) 

(92) 
(93) 
(94) 



These commutation relations forbid us from repeating the standard discussion of 1- 
dimensional shape invariance (see section |2|), and thus there is no automatic equivalence 
of the spectra for the Hamiltonian ( jS8"D and 



A' 



# (1) ("o) = E^(«o)4(«o) 



(95) 



i=l 
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However, if 

E A(a )4(ao) = E 4(ai)A(ai) + fl(c*i) , (96) 

the ground states of the two hamiltonians are related by 

^Woc^N. (97) 

Furthemore, the ground state of H^ has energy R(ai). 

Obviously, one can repeat the precedure, exactly as in the 1-dimensional case, and 
obtain the series of Hamiltonians 

n 

H<- n \a ) = # (0) K) + E R(<*k) ■ (98) 

fc=i 

The corresponding ground states are given by 

^ ft) (a )oc^ 0) K), (99) 

and have energies 

4 n) = E R(a k ) . (100) 

fc=i 

Unfortunately, there is no immediate relationship between these states and the spectrum 
of the original model. 

References 

[1] F. Cooper, A. Khare and U. Sukhatme. Supersymmetry and Quantum Mechanics. 
Phys. Rep., 251:267, 1995. 

[2] F. Calogero. Solution of a Three-Body Problem in One Dimension. J. Math. Phys., 
10:2191, 1969. 

[3] B. Sutherland. Exact Results for a Quantum Many-Body Problem in One Dimension. 
Phys. Rev., A4:2019, 1971. 

[4] M. A. Olshanetsky and A. M. Perelomov. Quantum Integrable Systems Related to 
Lie Algebras. Phys. Rep., 94:6, 1983. 

[5] A. Polychronakos. Exchange Operator Formalism for Integrable Systems of Particles. 
Phys. Rev. Lett, 69:703, 1992. 

[6] W. Ruhl and A. Turbiner. Exact Solvability of the Calogero and Sutherland Models. 
Mod. Phys. Let, A10:2213, 1993. 



22 



[7] A. Turbiner. Lie algebras and linear operators with invariant subspaces. In N. Kam- 
ran and J. Olver, editors, Lie Algebras, Cohomologies and New Findings in Quantum 
Mechanics, Contemporary Mathematics 160, page 263. AMS, 1994. 

A. Kostelecky, M. M. Nieto and D. R. Truax. Supersymmetry and the Relationship 
Between the Coulomb and Oscillator Problems in Arbitrary Dimensions. Phys. Rev., 
D32:2627, 1985. 

[9] D. Z. Freedman and P. F. Mende. An Exactly Solvable A^-Particle System in Super- 
symmetric Quantum Mechanics. Nucl. Phys., B344:317, 1990. 

[10] L. Brink, T. H. Hansson, S. Konstein and M. A. Vasiliev. The Calogero Model: 
Anyonic Representation, Fermionic Extension and Supersymmetry. Nucl. Phys., 
B401:591, 1993. 



[11] A. Das and S. A. Pernice. Higher dimensional susy quantum mechanics, [hep- 
I th/96l2l25|. 



[12] I. G. Macdonald. Symmetric Functions and Hall Polynomials. Clarendon Press, 
Oxford, 2nd edition, 1995. 

[13] Luc Lapointe and Luc Vinet. Exact Operator Solution of the Calogero-Sutherland 
Model. Comm. Math. Phys., 178:425, 1996. 

[14] F. Calogero. One-Dimensional Many-Body Problems with Pair Interactions Whose 
Exact Ground State Wave Function is of Product Type. Lett. Nuovo Cim., Ser. 2, 
13:507, 1975. 

[15] M. Abramowitz and I. A. Stegun. Handbook of Mathematical Functions with Formu- 
las, Graphs and Mathematical Tables. Dover, New York, 2nd edition, 1972. 



28 



